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ABSTRACT 

A comprehensive model atom was constructed for C I using the most up-to-date 
atomic data. We evaluated non-local thermodynamical equilibrium (NLTE) line for¬ 
mation for neutral carbon in classical ID models representing atmospheres of late-type 
stars, where carbon abundance varies from solar value down to [C/H] = —3. NLTE 
leads to stronger C I lines compared with their LTE strength and negative NLTE 
abundance corrections, Anlte- The deviations from LTE are large for the strong lines 
in the infrared (IR), with Anlte = —0.10 dex to —0.45 dex depending on stellar 
parameters, and they are minor for the weak lines in the visible spectral range, with 
|Anlte| ^ 0.03 dex. The NLTE abundance corrections were found to be dependent of 
the carbon abundance in the model. As the first application of the treated model atom, 
carbon NLTE abundances were determined for the Sun and eight late-type stars with 
well-determined stellar parameters that cover the —2.56 ^ [Fe/H] ^ —1.02 metallicity 
range. Consistent abundances from the visible and IR lines were found for the Sun 
and the most metal-rich star of our sample, when applying a scaling factor of S'!! = 
0.3 to the Drawinian rates of C-l-H collisions. Carbon abundances were also derived 
from the molecular CH lines and, for each star, they agree with that from the atomic 
C I lines. We present the NLTE abundance corrections for lines of C I in the grid of 
model atmospheres applicable to the carbon-enhanced (CEMP) stars. 

Key words: line: formation. Sun: abundances, stars: abundances, stars: late-type 


1 INTRODUCTION 

Carbon abundance studies are of high importance in vari¬ 
ous fields of astrophysics and astrochemistry. Carbon is an 
essential element for the beginnings of life on the Earth. It 
acts as a primary catalyst for H-burning via the CNO cycle, 
and it contributes significantly to the stellar interior and at¬ 
mosphere opacity. Carbon plays an important role in dust 
formation processes in the interstellar medium. 

Carbon is one of the representatives of the ’heavy’ 
{Z ^ 6) elements, which are of stellar nucleosynthesis ori¬ 
gin. Accurate stellar carbon abundances may improve our 
understanding the chemical evolution history of our Galaxy 
and other galaxies. In the course of normal stellar evolution, 
carbon is essentially all produced by He burning through th e 
3a nuclear reaction: 3"^He —^ ^^C -I- 7 llBurbidge et al.lll957ll . 
In the modern galacto-chemical evolution models the major 
sources of carbon are supernovae of type H (SNell), hyper¬ 
novae (HNe), supernovae of type la (SNela) , and asymptotic 
giant branch (AGB) stars (see, for example lKobavashi et aH 


[ 2 ^. A relative contribution of different sources to the 
galactic carbon varied with time. 

In the early Galaxy, the carbon production was domi¬ 
nated by SNell and HNe. Indeed, progenitors of SNell and 
HNe are massive stars, with a mass of larger than 8 and 
20 solar mass (Mq), respectively, and a short lifetime of few 
10®-10^ years. One also considers the scenarios with rotating 
massive stars. The delay time for nucleosynthesis in SNela 
ranges between 0.3 Gyr and 3 Gyr, for a w ide variety of 
hypothesis on th e prog enitors (lGreggidl2005l '). According to 
iKobavashi et al.1 (l201lll . an onset of the carbon production 
by the AGB stars refers to the time, when the galactic iron 
abundance grew to [Fe/H] ~ —1.5. 

With the SNell -|- HNe model, IKobavashi et al.] ll201ll ~) 
predicted the close-to-solar G/Fe abundance ratio at low 
metallicities, —4 ^ [Fe/H] Si —3, but the [C/Fe] ratio be¬ 
comes as large as 0.9, when including the yields of rotating 
massive stars. At [Fe/H] ~ —1, [G/Fe] reaches 0.13, if all the 
sources, AGB -I- SNela -I- SNell -I- HNe, are included, and 
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[C/Fe] decreases down to —0.19, when excluding the AGB 
yields. 

A role of different carbon sources at different metallic- 
ities can be investigated from confronting the observations 
with the chemical evolution models. Studies in the litera¬ 
ture demonstrate a variety of the observed trends of [C/Fe] 
Fe/Hj. For the s ample of the —3.2 < [Fe/H] < —0.7 
Akerman et al.l ll2004h obtained enhanced C abun- 


versus 
dwarfs 

dance relative to Fe, with [C/Fe] reducing towards higher 
metallicity from 0.45 to 0.25. They analysed the near in¬ 
frared (IR) C I lines unde r the local thermody namic equi¬ 
librium (LTE) assumption. iFabbian et ^ ll2006ll revised the 
C abundances of that stellar sample based on the non-local 
thermodynamic equilibrium (NLTE) line formation for C I 
and obtained the [C/Fe] ratio to be, on average, close to 
the solar one. iReddv et al.l (l2006l l studied the C I lines in 
the visible spectral range in the sample of the thick disk 
and thin disk F-C dwarfs in the —1.2 < [Fe/H] < -1-0.2 
range. From the LTE analysis they obtained a supersolar 
[C/Fe] ratio of about 0.4 for the most metal-poor ([Fe/H] 

< —0.4) thick disk stars and a decline down to the subsolar 
values at higher metallicity. In the overlapping metallicity 
range, the thin disk stars reveal, on average, lower [C/Fe] 
ratios than the thick dis k stars. In contrast, the study by 
iBensbv fc Feltzin3 ll2006h based on the [C I] line at 8727 A 
shows that [C/Fe] versus [Fe/H] trends for the thin and thick 
discs are totally merged and flat, with [C/Fe] ~ 0.1, for the 
subsolar metallicities, down to [Fe/H] = —0.8. 

From observati ons of the CH bands in 83 subdwarfs 
ICarbon et al.l ll 19871 1 found [C/Fe] to be essentially constant 
and close to the solar value over the range —2.5 < [Fe/H] 

< —0.7. However, they noted an upturn in the [C/Fe] values 
at [Fe/H] < —2. A remarkably fiat [C/Fe] versus [Fe/H] rela, - 
tion, with [C/Fe] = -1-0.18, was found bv ISpite et al.l (l2005l l 
from observations of the CH bands for the sample of ’un¬ 
mixed’ giant stars in the range —4 < [Fe/H] < —2.5. 

The situation with understanding sources of the 
carbon production in the early Calaxy is complicated 
by a discovery of the Carbon-Enhanced Metal-Poor 
(CEMP) stars that show very high enhancements in 
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ments are inborn, or their atmospheres were polluted, most 
likely by accretion from an ACB b inary companion. Recent 
study of IStarkenburg et ^ ll2014ll conhrmed a binarity of 
the CEMP-s stars that are additionally enhanced in barium. 
But new radial velocity data for the 15 CEMP-no stars, not 
enhanced in barium, were found to be inconsistent with the 
binary properties of the CEMP-s class, thereby strongly in¬ 
dicating a different physical origin of their carbon enhance¬ 
ments. 

Stellar carbon abundance determinations rely on vari¬ 
ous spectroscopic indicators. These are allowed and forbid¬ 
den lines of C I and lines of the molecular species CH, C 2 , 
and CO. Suitable lines of C I are located in the visible (4300 
- 7900 A) and near-IR (7900 - 20000 A) spectral regions. All 
the allowed lines have close together excitation energies of 
the lower level, Eexc, but different oscillator strengths, with 


smaller values for the visible than the near-IR lines. As a re¬ 
sult, the visible C I lines are much weaker than the near-IR 
ones. 

The [C I] forbidden line at 8727 A can reliably be mea¬ 
sured for the close-to-solar metallicity stars. For stars, not 
enhanced in carbon, the C I visible lines can be used down 
to [Fe/H] = —1.5 and the IR lines down to [Fe/H] = —2.5. 
At the lower metallicity, [Fe/H] < —2.5, stellar carbon abun¬ 
dances are determined presumably from the CH molecular 
lines. 

An use of different abundance indicators with classi¬ 
cal ID model atmospheres can produce systematic shifts 
and false metallicity tr ends. This is exactly the case in 
the lTomkin et ahl (Il992fl study. They obtained, on average, 
0.4 dex higher abundance from the atomic C I compared 
with the molecular CH lines for the sample of metal-poor 
halo dwarfs. When dealing with a wide metallicity range, it 
is impossible to apply a common carbon abundance indica¬ 
tor. One needs, therefore, to investigate the possible sources 
of abundance discrepancy between the atomic and molecu¬ 
lar lines. These can be (i) employing the homogeneous and 
plane-parallel (ID) model atmospheres, (ii) inadequate line- 
formation treatment based on the LTE assumption, (iii) and 
the uncertainties in stellar parameters. 

First determinations of the carbon abundance based 
on a three-dimensional (3D), time-depende nt and hydro- 
dynam ical model atmosphere were made bv lAsplund et ahl 
(l2005ll for the Sun. For the atomic lines, it was shown that 
the abundance differenc es between the 3D and ID MARCS 
llCustafsson et aP Il975ll models are small and amount to 
-1-0.01 dex, -fO.Ol dex, and —0.02 dex, on average, for the 
atomic C I, mole cular CH elect r onic and C 2 electronic 
lines, respectively. lAsplund et ahl (l2009ll updated the cal¬ 
culations and obtained consistent abundances from the dif¬ 
ferent species, with the mean log ec = 8.43 ± 0.05 for the 
3D mo del and log ep = 8.42 ± 0.05 for the ID MARCS 
model. ICaffau et ahl (l2010fl computed positive 3D-ID abun¬ 
dance corrections for the solar C I lines, up to 0.1 dex for the 
stronger near-IR lines and ^0.03 dex for the fainter visible 
lines. 

For the stellar atmosphere parameters beyond the so¬ 
lar ones the 3D ca lculations were only performed for the 
Cl fictitious lines (jPobrovolskaj [20131 . iDobrovolskas et al.1 
I2OI3II . For the Eexc = 6 eV lines at A = 8500 A in 
the Teff/log 5 /[M/H] = 5930/4.0/0, 5850/4.0/-1, and 
5860/4.0/—2 models the 3D-1D abundance corrections do 
not e xceed few 100-th in a logarithmic scale (jPobrovolskasI 
l2013l. Fig. 2.19). 

An influence of the 3D effects on abundance deter- 
mination from th e m olecular CH l ines w as evaluated by 
ICollet et al.1 ll2007ll andiHav ek et ah l laoTil) for cool giants of 
various metallicity. ICollet et al. 1 20o5'^howed that the 3D- 
ID abundance corrections for the CH lines at A ~ 4360 A 
are overall negative and they amount to —0.15 dex at so¬ 
lar metallicity and —0.8 dex at [M/H] = — 3 in the models 
with Teff = 5050 K and log q = 2.2. iHavek et al.1 (1201 ill 
updated the ICollet et al.1 (l2007l i calculations and obtained 
smaller 3D effects, namely 3D—ID = 0.0, —0.18 dex, and 
—0.35 dex in the giant (log g = 2.2) models with [M/H] = 
0 (Teff = 5060 K), [M/H] = -2 (5050 K), and [M/H] = -3 
(5100 K), respectively. 

Thus, the discrepancy between the C I and CH-based 
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abundances found by iTomkin et alJ lll992ll in the classical 
ID analysis is expected to be larger in 3D. 

The NLTE m eth ods for C I were developed 
by IVe rnazza et al.| l|l981^, IShchukina fc Shcherbinal 
90D, ' 


(Il990li. Stuerenburg fc Holweger (ll99Clll. Ta,kedal ( 1992ll . 


Prz^MUa_e^^alJ_^ 200lj)' iFabbian et al. (l200y) , and 


Lvubimkov et al.l (l2015ri . IStuerenburg fc Holweged (Il990l l 


showed that NLTE removes a discrepancy between visible 
and near-IR lines of C I in the Sun, while the difference 
in LTE abundances between individual lines can be up 
to 0.37 dex, for exampl e between C I 9111 A and 5380 A. 
iTakeda fc Hondal ll2005l i and IFabbian et al.l ll2006li inves¬ 
tigated the deviations from LTE for C I in wide stellar 
parameter range and found that the NLTE corrections are 
overall negative and grow in absolute value towards lower 
carbon abundance. However, these two studies provide dif¬ 
ferent results for hot very metal-poor models, for example 
the difference in NLTE correction for C I 9094 A amounts to 
0.15 dex in the 6000/4.0/—2 model. IFabbian et al.l (l2006l l 
showed that an adequate line formation modelling for C I is 
important for stellar abundance determinations in Galactic 
chemical evolution studies, where the stellar sample covers 
a range o f metallicities of more t han 2 dex. Indeed, the LTE 
results of lAkerman et aP (l2004l i revealed an Q-element like 
behavior of carbon in the —3 Si [Fe/H] ^ —0.7 metallicity 
range, while c l ose-to -solar C/Fe ratios were obtained by 
IFabbian et al.l ll2006tl after correcting these data for the 
NLTE effects. The Lvubimkov et all (l2015h studied the 
NLTE effects for C I in the galactic A-, F- and G-type 
supergiants and bright giants. 

This study aims to construct a comprehensive model 
atom for C I using the most up-to-date atomic data avail¬ 
able so far and to determine the carbon abundance of the 
Sun and selected late-type stars with well-determined stel¬ 
lar parameters based on the NLTE line-formation for G I in 
the classical ID model atmospheres and high-resolution ob¬ 
served spectra. Our small sample covers the —2.58 ^ [Fe/H] 
Si 0 range. We also derive the carbon abundance from the 
molecular GH and C 2 lines and investigate the differences 
between the atomic and molecular lines, which can serve to 
characterise empirically how neglecting the 3D effects influ¬ 
ences the derived abundances. 

This paper is organized as follows. The model atom of 
carbon, the adopted atomic data, and the departures from 
LTE for C I depending on stellar parameters are described 
in Sect.[2l Analysis of the atomic G I and molecular CH 
and G 2 lines in the Sun is presented in Sect.O In Sect.|4l 
we determine the C abundance of the selected stars. We 
summarise our conclusions in Sect .(5] 


0-2), 2s^2p nf (n = 4-8), and 2s2p® electronic configurations 
and the 2s2p® quintet term. In addition, the 2s^2p nl (nl = 
5g, 6g, 6h) levels were taken from the Kurucz’s databasfl 
To provide close collisional coupling of G I to the contin¬ 
uum electron reservoir, the energy separation of the highest 
G I levels from the ionization limit must be smaller than 
the mean kinetic energy of electrons, i.e., 0.5 eV for atmo¬ 
spheres of solar temperature. In our model atom the energy 
gap between the uppermost levels of G I and ionization limit 
is equal to 0.13 eV. Fine structure splitting was included ev¬ 
erywhere, up to n = 8. All the states with n = 9 and n = 10 
have close energies, and their populations must be in equilib¬ 
rium to each other. Therefore, the levels of common parity 
were combined into the single superlevel. The Grotrian term 
diagram for our model atom is shown in Fig.[T] 

Radiative data. Our model atom includes 1524 allowed 
bound-bound (6 — b) transitions. Their transit ion probabil¬ 
ities were taken from the NIST and VALD llKupka et al.l 
I 1999 II databases, wher e available, and th e Opacity Project 
database TOPbas^l (iGunto et al.l Il993l . ILuo et Si Il989l . 
iHibbert et al.l[l993fl . Photo-ionization cross-sections for lev¬ 
els with n ^ 8, I ^ 3 were from TOPbase, and we adopted 
the hydrogen-like cross-sections for the two superlevels. 
Collisional data. All levels in our model atom are cou¬ 
pled via collisional excitation and ionisation by electrons 
and by neutral hydrogen atoms. Detailed electron-impact 
excitation cross-sections calculated in the close-coupling ap¬ 
proximation using the R-matrix meth o d are available for 
more than 400 transitions from iRei ^ (I 1994 I L For the re- 
maini ng transitions we use the impact parameter method 
(IPM. [Seaton|[l962bll for the allowed transitions and assume 
that the effective collision strength Hij = 1 for the forbid¬ 
den transitions. Accurate data on inelastic collisions of car¬ 
bon with neutral hydrogen atoms remain still undefined. W e 
employed, therefore , the Drawin’s formula (|DrawirJ ll968l L 
as implemented by Isteenbock fc Holwegeil (ll984l L The ef¬ 
ficiency of G-l-H collisions is treated as a free parameter 
in our attempt to achieve consistent element abundances 
derived from the visible and near-IR lines in the Sun and 
selected stars. For each object, the calculations were per¬ 
formed with a scaling factor of Sn = 0, 0.1, 0.3, and 1. Hy¬ 
drogen collisions for the forbidden transitions were ignored. 
Ionisation b y electr o nic col lisions was everywhere treated 
thro ugh the SeatorJ (ll962a|) clas sical path approximation. 
The lSteenbock fc Holweger! (Il984l ') formula, with Sh was ap¬ 
plied to calculate the C-l-H ionisation rates. _ 

_ Our model atom is similar to that of IFabbian et al.l 

1 I 2 OO 6 II except applying in this stu dy d e tailed electron- 
impact excitation cross-sections from iRei dl (ll994l) . 


2 NLTE LINE FORMATION FOR C I 

2.1 Model atom and atomic data 

Energy levels. Model atom includes 208 energy levels of 
G I up to n = 10, Z = 4, nine lowest levels of C H, and the 
ground state of G HI . Most energy levels w ere taken from 
the NIStQ database (iRalchenko et al]l2008l ). They belong 
to singlet and triplet terms of the 2s^2p nl (n = 2-10, I = 


2.2 Method of calculations 


To solve the radiative transfer and sta tistical equilibrium 
equa tions, we used the code DETAIL ilButler fc Giddingj 
Il985ll based on the a ccelerated A-iteration method 
IRvbicki fc Hummer! Il991^. The opacity package was im¬ 


proved as described in iMashonkina et al 


y pacK a 

llml). 


The depar¬ 


ture coefficients, 6i = riNLTE / were then used to com- 


^ http://physics.nist.gov/PhysRefData/ 


^ http://cfaku5.cfa.harvard.edu/atoms.html 
^ http//legacy.gsfc.nasa.gov/topbase 
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Figure 1. Grotrian term diagram for neutral carbon. The solid lines indicate the seven transitions, where the investigated 
spectral lines arise. The nf, ng, nh energy levels, where the Russell-Saunders coupling is broken, are shown on the right 
column. 


pute synthetic line profiles via t he SIU (Spectrum Investiga¬ 
tion Utility) code (^eeti[l993). Here, 51-nlte and wlte are 
the statistical equilibrium and thermal (Saha-Boltzmann) 
number densities, respectively. 

The calculations were performed with using plane- 
parallel, homogeneous (ID) MARCS model atmospheres 
llGustafsson et al]l2008l '). 


2.3 Departures from LTE depending on carbon 
abundance and surface gravity 

Figure[2] shows the departure coefficients for the selected 
levels in the solar 5777/4.44/0 model atmosphere. Neutral 
carbon C I dominates the element number density over all 
atmospheric depths. Thus, no process seems to affect the 
C I ground-state and low-excitation level populations sig¬ 
nificantly, and they keep their thermodynamic equilibrium 
values. The levels 3s^P° and 3s®P° are overpopulated out¬ 
ward logT = —0.2 due to the ultraviolet radiative pumping 
in the transitions 2p^ - 3s^P° (2479 A), 2 p^ - 3s^P° 

(1931 A), and 2p^ ®P - 3s®P° (1657 A). Levels with Eexc 
> 8 eV are underpopulated outward logr = —0.2 due to 
spontaneous transitions to the lower levels. 


The NLTE effects for a given spectral line can be under¬ 
stood from analysis of the departure coefficients at the line 
formation depths. The C I lines used in abundance analysis 
are listed in Table[T] Here, we consider two lines with sim¬ 
ilar Eexc, but different ^/-values. The C I 5380 A (3s^P° - 
4p^P) line is weak, and it forms in the deep layers, around 
logrsooo = 0, where the departures from LTE are small. In 
contrast, C I 9111 A (3s^P° - 3p®P) is strong, and its core 
forms around logrsooo = — 1, where the departure coefficient 
of the lower level, fe(3s®P°) is larger than unity and fe(3s®P°) 
> 6(3p^P). NLTE leads to strengthened lines of C I and neg¬ 
ative abundance corrections, Anlte = logCNLTB - logeuTE, 
with the stronger NLTE effects for the stronger lines. 

To evaluate the NLTE effects depending on the car¬ 
bon abundance, we calculated the NLTE abundance correc¬ 
tions with the four model atmospheres having common Tefi 
= 6000 K and logg = 4.0, but different metallicities, with 
[M/H] = 0, —1, —2, and —3. For each model the NLTE 
calculations were performed with the two different carbon 
abundances, [C/Fe] = 0 and [C/Fe] = 1. The Sh = 0.3 and 
a microturbulence velocity of = 1 Lms”*^ were employed. 
The obtained results are displayed in Fig. |3] for C I 9405 A 
and 9111 A. 
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Figure 2. Departure coefficients for the C I levels and the C II 
ground state as a function of logrsooo in the solar MARCS 
5777/4.44/0.00 model atmosphere. The vertical lines mark the 
line-center formation depths for C I 9111 A and C I 5380 A. 




Figure 3. NLTE abundance corrections versus line equivalent 
width (EW) for C I 9405 A (solid curve) and 9111 A (dashed 
curve) in the model atmospheres with Tgff = 6000 K, logg = 
4.0, and different metallicities of [M/H] = 0, —1, —2, and —3 
(circles, squares, up-triangles, and down-triangles respectively). 
The NLTE calculations were performed with [C/Fe] = 0 (left 
panel) and [C/Fe] = 1 (right panel). Everywhere, Sjj = 0-3 and 
= 1 kms“t. 


The NLTE effects for C I depend on the element abun¬ 
dance. For [C/Fe] = 0 the NLTE abundance corrections 
reduce, in absolute value, towards lower equivalent width 
(element abundance) throughout the metallicity range from 
[M/H] = 0 to —3. To understand this, we consider the depar¬ 
ture coefficients in Fig.[4l On the one hand, the total number 
of free electrons, acting as the source of thermalisation, de¬ 
creases with decreasing metallicity and the departures from 
LTE are amplified. Indeed, at a common optical depth the 
departure coefficients deviate from unity larger in the low 
than in the solar metallicity model. On the other hand, with 
metallicity (carbon abundance) decreasing, the line forma¬ 
tion region moves to deeper layers, where the NLTE effects 
wane. Of the two competing effects the latter prevails, and 
the NLTE effects reduce toward lower metallicity. 

A different behavior was found in the [C/Fe] = 1 case, 
namely the NLTE effects grow, when moving from [M/H] 
= 0 to —1, and then reduce towards lower EW. The C I 
9405 A and 9111 A lines are very strong in the [M/H] = 0 
model, with EW ~ 500 mA and 330 mA, respectively, which 
are mostly contributed from the broad line wings. The line 
wings form in deep atmospheric layers, where the departures 
from LTE are small. When moving to [M/H] = — 1, a con¬ 
tribution of the line wings decreases, and the NLTE effects 
grow. In the [M/H] = —2 and —3 models the NLTE effects 
decrease due to shifting the line formation regions in the 
deeper layers. 

We compared our calcula tions with the NLTE re - 
sults of iTakeda fc Hond^ (|2005l ) and iFabbian et ahl l|200(i ). 
To be as clo s e as possible to the NLTE method of 
IFabbian et_ajj_l[2 n0(jl^ we em ployed the same collisional 


Fabbian et al.l ll200fil ) and Sh = 1. Figure[5] shows 


recipe as 

the NLTE abundance corrections depending on metallicity 
for the representative line at 9094 A in the model atmo¬ 
spheres with common Teff = 6000 K, logg = 4.0, [C/Fe] = 
0.4, and = 1 kms“^. All the three studies give consistent 
results for the [M/H] = 0 and —1 models, where the NLTE 
effects ar e the strongest. A t [M/H ] = —2 our data agree 
well with iTakeda &: Hondal (l2005l) . while a discrepancy o f 
0.08 dex in Anltb was obtained wiH]jFabbhm_et . 

The NLTE correction computed bv iFabbian et all ( 20061 ) for 
the [M/H] = —3 model is as large as that for [M/H] = —2, 
although the C I 9094 A line is much weaker in the [M/H] 
= -3 (EW = 8 mA ) than [M/H] = -2 (EW = 50.2 mA ) 
model. In our calculations Anlte = —0.07 dex. It is worth 
noting that spectral lines of EW = 8 mA are too weak for ac¬ 
curate measurements, and, thus, the theoretical predictions 
cannot be checked with the observations. 


3 SOLAR CARBON ABUNDANCE 
3.1 Atomic C I lines 

As a first application of the treated model atom, we de¬ 
rived the solar carbon abundance from lines of C I. For com¬ 
parison, the element abundance was also determined from 
the molecular CH lines. The solar fl ux observations wer e 
taken from the Kitt Peak Solar Atlas llKurucz et al.l[l983 ). 
We used the MARCS model atmosphere 5777/4.44/0 and a 
depth-independent microturbulence of 0.9kms“^. The ele¬ 
ment abundance was determined from line profile fitting. As 
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Figure 4. Departure coefficients, 6, for the C I levels and the ground state of C II as a function of log r5ooo in th® four model 
atmospheres with different metal abundance: [M/H] = 0, —1, —2, and —3. Everywhere, = 6000 K, logg = 4.0, [C/Fe] = 1, = 

1 kms“^. In each panel the two vertical lines indicate the formation region for C I 9405 A. 



Figure 5. NLTE abundance corrections d epending on m e tallic- 
ity for C I 9094 A from this stu dy (cir cles), [Fabbian et al.l ll2006l . 
squares), and iTakeda &; Hondal ll2005l . rhombi). Everywhere, T'efr 
=6000 K, log g = 4.0, [C/Fe] = 0.4, = Ikms-i, 5 h = 1. 


a rule, the uncertainty in fitting the observed proHle is less 
than 0.02 dex for weak lines and 0.03 dex for strong lines. 
Our synthetic flux profiles were convolved with a profile that 
combines a rotational broadening of 1.8 kms“^ and broad¬ 


ening by macroturbulence with a radial-tangential profile. 
The most probable macroturbulence velocity Vmac was var¬ 
ied between 2kms“^ and 4kms“^ for different lines of C I 
and CH. Quality of the hts is illustrated in Fig. [6] for two 
lines of C I. The C I 5380 A profile beyond 5380.45 A is 
contributed from an unknown blend. Weak telluric line at 
9111.95 A affects C 1 9111 A, however, the difference be¬ 
tween observed and calculated NLTE spectra, (O - C), does 
not exceed 0.4 % for the remaining prohle. 

For lines listed in Table[T] we determined the element 
abundance under various line-formation assumptions. In 
LTE, the abundance difference between the visible and near- 
IR lines of C I, Aloge(vis - IR), amounts to —0.21 dex. In 
line with the previous studies, we Hnd that the [C I] 8727 A 
forbidden line does not suffer from the NLTE effects, be¬ 
cause it arises in the transition between the metastable lev¬ 
els that keep the LTE populations. The NLTE corrections 
are small for the visible lines, independent of the applied col- 
lisional recipes, with |Anlte| ^ 0.03 dex. For the IR lines, 
the departures from LTE are sensitive to a variation in col- 
lisional rates. For example, for different lines Anltb ranges 
between —0.14 dex and —0.45 dex, when Sh = 0, and be¬ 
tween —0.08 dex and —0.32 dex, when 5 h = 1. Consistent 
within 0.03 dex and 0.02 dex NLTE abundances from the 
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Table 1. Lines of C I and solar carbon abundance from atomic lines 


A, A 

Transition 

Mult. 

log gf 

-^exc 

LTE 


NLTE, 

, S'h = 









0.0 

0.1 

0.3 

1.0 

8727.139 

2p2 M 

- 2p2 Is 


-8.165 

1.26 

8.45 

8.45 

8.45 

8.45 

8.45 

4932.049 

3sip? 

- 4p®So 

13 

-1.658 

7.69 

8.48 

8.45 

8.45 

8.45 

8.45 

5052.167 

Os^P? 

- 4p^D2 

12 

-1.303 

7.69 

8.50 

8.47 

8.47 

8.47 

8.47 

5380.337 

Os^P? 

- 4piPi 

11 

-1.616 

7.69 

8.46 

8.43 

8.43 

8.43 

8.43 

6587.610 

3piPi 

- 4dipS 

22 

-1.003 

8.54 

8.45 

8.43 

8.43 

8.43 

8.43 

Mean (vis) 






8.47 

8.45 

8.45 

8.45 

8.45 

a 






0.02 

0.02 

0.02 

0.02 

0.02 

8335.148 

Os^P? 

- 3piSo 

10 

-0.437 

7.69 

8.54 

8.40 

8.40 

8.43 

8.46 

9078.288 

3s® P? 

- 3p®Pi 

3 

-0.581 

7.48 

8.65 

8.41 

8.44 

8.46 

8.51 

9094.834 

3s®P^ 

- 3p®P2 

3 

0.151 

7.49 

8.72 

8.34 

8.37 

8.40 

8.46 

9111.809 

3s® P^ 

- 3p®Pi 

3 

-0.297 

7.49 

8.67 

8.34 

8.37 

8.40 

8.46 

9405.730 

3s®P)’ 

- 3p^D2 

9 

0.286 

7.69 

8.79 

8.34 

8.37 

8.40 

8.47 

9658.431 

3s®P5 

- 3p®Si 

2 

-0.280 

7.49 

8.71 

8.37 

8.40 

8.43 

8.48 

9061.433 

3s® P? 

- 3p®P2 

3 

-0.347 

7.48 

- 

- 

- 

- 

- 

9062.492 

3s® PS 

- 3p®Pi 

3 

-0.455 

7.48 

- 

- 

- 

- 

- 

Mean (IR) 






8.68 

8.37 

8.39 

8.42 

8.47 

a 






0.08 

0.03 

0.03 

0.02 

0.02 

Mean 






8.59 

8.41 

8.42 

8.43 

8.46 

CT 






0.12 

0.05 

0.04 

0.03 

0.02 


Mult.: the multiplet numbers accordingly to^oore Q,972). 


visible and near-IR lines were obtained for Sh = 0.3 and 1, 
respectively. We derived from the atomic lines a solar carbon 
abundance of logec = 8.43±0.03 (S'h=0.3). 


3.2 Molecular C 2 and CH lines 


C 2 electronic lines. We selected the 12 least blended 
lines from the C 2 Swan band in the 4992-5150 A range. 
They are listed in Table[2] along with the lower level exci- 
tation potent i als an d the transition probabilities taken from 
I Brooke et al.l (l2013ll . 

CH electronic line s. We used a dissoci ation energy of 
Do(C 2 ) = 6.297 eV from IUrdahl et all lll99ll ~). Our analysis 
is based on the lines from the (0, 0) and (1, 1) bands of CH 
A-X around 4300 A. They are listed in TableO The lower 
level excitation potentials and the tr ansition probabilities 
were taken from iMasseron e t al. (l20 14h. The dissociati on en¬ 
ergy is Do(CH) = 3.465 eV I Huber fc Herzberglll979l i. Our 
test calcul ations showed that usi ng the molecular param¬ 
eters from I Jorgensen et al.l (Il996li results in higher abun¬ 
dances, by 0.00 to 0.04 dex for different CH lines, compared 
with t he corresponding values based on the lMasseron et al.l 
(l2014h data. 


Abundances from the individual molecular C 2 and CH 
lines are presented in Tables[2] and [3] It can be seen that 
lines of C 2 give consistent within 0.04 dex abundances, and 
the mean abundance, log ec = 8.46±0.02, agrees well with 
the NLTE(S'h = 0.3) abundance from the atomic C I lines. 
Abundances from different lines of CH are consistent within 
0.05 dex, and the mean abundance, log ec = 8.39T0.02, 
agrees with that from the atomic lines. 


3.3 Comparison with previous studies 


When applying a common Sh of 0, we find the obtained solar 
mean abundance from lines of C I (8.41T0.05) to be consis- 
tent within th e erro r bars with the ID and 3D results of 
lAsplund et al.l (l2009l . hereafter, AGSS09), namely log ec = 
8.39±0.04 (ID) and 8.42±0.05 (3 D), and the 3D abundance 
log ec = 8.446T0.121 calculated bv ICaffau et al.l (l2010l . here¬ 
after, CLBIO). Abundance from the forbidden [C I] 8727 A 
line, log ec = 8.45, is 0.07 dex and 0.10 dex higher compared 
with the corresponding values from AGSS09 and CLBIO. 
In part, this is due to employing in this study a 0.03 dex 
lower g f-value . Indee d, we relied on log gf = —8.165 from 
I Froese F ischeil (l2006ll. while AGSS 09 and CLBIO adopted 
log gf = —8.136 bv lHibbert et al.l (ll993l L Another source 
of the discrepancy is a diff erent treatmen t of faint blend¬ 
ing lines. According to the iKuruc i (l2007ll data, two lines 
of Fe I contribute to the 8727 A blend. These are Fe I 
8727.10 A (Sexc = 5.587 eV, log gf = -5.9 24) and Fe I 
8727.1 3 A (Sexc = 4.186 eV, log gf = -4.262). [CafeueTaD 
ll2010lf treated Fe I 87 27.132 A , using the higher log gf = 
—3.93 from the older iKurucd (Il993l l calculations. A dif¬ 
ference of —0.33 dex in gj-value of the blending Fe I line 
makes a -1-0.02 dex difference in the carbon abundance de¬ 
rived from [C I] 8727 A. To understand a source of the 
remaining dis crepancy, we inspec ted the forbidden line in 
the solar flux (iK urucz et al.l ll984l l and disk-center intensity 
llBrault fc Testerman 1972l l spectra and obtained fully con¬ 
sistent abundances, with logec = 8.45. It is worth noting, 
we measured EW([C I] 8727) = 5.5 mA in the disk-center in¬ 
tensity spectrum, while AGSS09 and CLBIO reported lower 
values of 5.3 mA and 5.1 mA, respectively. 

Our ID results for the molecular lines agree well with 
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Table 2. Carbon abundances derived from solar lines of C 2 


A, A 

log(3/) 

E 

-^exc 

loge(C) 

4992.2750 

0.288 

0.802 

8.44 

4992.3035 

0.281 

0.802 

5033.7792 

0.191 

0.584 

8.47 

5143.3240 

-0.411 

0.102 

8.44 

5144.9149 

-0.447 

0.097 

8.46 

5145.2255 

-0.485 

0.098 

8.46 

5150.5448 

0.026 

0.328 


5150.5542 

0.038 

0.328 

8.46 

5150.6461 

0.013 

0.328 

5150.6737 

0.638 

0.084 


5052.6161 

0.153 

0.489 

8.45 

5052.6254 

0.162 

0.489 

5086.3897 

0.031 

0.328 

8.48 

5103.7231 

-0.023 

0.250 

8.47 

5103.7710 

-0.038 

0.250 

5109.0921 

-0.053 

0.227 


5109.1490 

-0.068 

0.228 

8.47 

5109.3030 

-0.084 

0.228 


5135.5509 

0.137 

0.471 


5135.5818 

0.127 

0.472 

8.45 

5135.6825 

0.117 

0.472 


5073.4490 

0.090 

0.388 


5073.4513 

0.101 

0.388 

8.45 

5073.5815 

0.080 

0.388 


Mean 



8.46±0.02 


Table 3. 

Carbon abundances derived from solar 

A, A 

log(9/) 

Eexc[eV] 

logec 

4218.710 

-1.315 

0.413 

8.37 

4218.734 

-1.337 

0.413 

4248.729 

-1.467 

0.191 


4248.937 

-1.431 

0.191 

8.38 

4248.952 

-3.256 

0.191 


4253.000 

-1.506 

0.523 

8.41 

4253.206 

-1.471 

0.523 

4255.248 

-1.461 

0.157 

8.41 

4255.248 

-3.210 

0.157 

4263.969 

-1.575 

0.459 

8.36 

4274.133 

-3.025 

0.074 

8.38 

4274.186 

-1.563 

0.074 

4356.355 

-1.846 

0.157 

8.40 

4356.371 

-1.455 

1.109 

4356.594 

-1.793 

0.157 

8.39 

Mean 



8.39±0.02 



W avelength( Angstrom) 



W avelength( Angstrom) 


Figure 6. Best NLTE fits (continuous curve) of the solar C I 
5380 A and 9111 A lines (asterisk). For each line, the LTE pro¬ 
file (dashed curve) was computed with the carbon abundance ob¬ 
tained from the NLTE analysis. The differences between observed 
and calculated NLTE spectra, (O - C), are shown in the upper 
parts of the panels. 


the ID and 3D data of AGSS09, log ec = 8.40±0.03 (CH) 
and 8.46±0.03 (C 2 ) in ID and log ec = 8.43±0.03 (CH) and 
8.46±0.03 (C 2 ) in 3D. 


4 DETERMINATION OF CARBON 

ABUNDANCES OF THE SELECTED STARS 

4.1 Stellar sample, observations, and stellar 
parameters 

We selected nine metal-poor stars in the —2.56 < [Fe/H] 
^ —1.02 metallicity range, for which the high-quality ob¬ 
served spectra and the well determined stellar parameters 
are available in the literature (Tabled . We describe briefly 
the sources of stellar effective temperatures and surface grav¬ 
ities. 

Effective te mperatures. For HD 1225 63 we adopted 
Teff = 4600 K from liviashonkina et al.l (1201 ll ). and it is con¬ 
sistent with the recent value based on the interferometric ob¬ 
servations of ICreevev et al.l (l2012l) . The infrared flux method 
(IRFM) te mperatures were adopted for HD 84937 (as recom¬ 
mended by iMashonkina et al.l[20lll ). HD 59374, HD 94028, 
HD 103095, HD 140283, BD-4°3208, and BD-b66°268, as 
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Table 4. Atmospheric parameters of the selected stars and 
sources of the data. 


Star 

Teff 

[K] 

log 9 
(CGS) 

[Fe/H] 

[kms“^ 

Ref. 

1 


Reference stars 



HD 29907 

5500 

4.64 

-1.55 

0.6 

1 

HD 59374 

5850 

4.38 

-1.02 

1.2 

3 

HD 84937 

6350 

4.09 

-2.08 

1.7 

2 

HD 94028 

5970 

4.33 

-1.50 

1.4 

3 

HD 103095 

5130 

4.66 

-1.26 

0.9 

3 

HD 122563 

4600 

1.60 

-2.56 

2.0 

2 

HD 140283 

5780 

3.70 

-2.38 

1.6 

3 

BD-4°3208 

6390 

4.08 

-2.20 

1.3 

3 

BD+66°268 

5300 

4.72 

-2.06 

0.7 

3 

Notes. References: 

IMashonkina et alJ i 

2^; (2) 

IMashonkina et alJ il201llk f3) ISitnova et alJ ||2015 



recommended by 

Sitnova et al.l (120151) based on the data 

of Alonso et al.l ( 

1996), 

Gonzalez Hernandez & Bonifacio 

I2009f). ICasagrande et al. 

(2010l). and ICasaerande et al. 

120111). For HD 29907 its effective temperature was deter- 


mine d from the Balmer line wing fits (|Mashonkina et ^ 
l2003fl . 


Surface gravities. For HD 29907, HD 84937, 
HD 122563, and HD 140283 we adopted the Hip- 
parcos parallax-based surfac e gra vities calcula t ed b y 
iMashonkina et aP ll2003l . [20l3), and ISitnova et al.l (l2015ll . 
For the remaining five stars, log g were derived by 
ISitnova et alJ (12013 1 from the NLTE analysis of the Fe I 
and Fe H lines. 

Spectral observations for four program stars were car¬ 
ried out with the UVES spectrograph at the 8-m VLT tele¬ 
scope of the European Southern Observatory (ESO) and for 
five stars with the Shane 3-m telescope of the Lick observa¬ 
tory using the Hamilton echelle-spectrograph. Characteris¬ 
tics of the observed spectra are summarized in TableO All 
spectra were obtained with a high spectral resolving power, 
R, and high S/N ratio (more than 100). For HD 84937, 
HD 122563, and HD 140283 their observe d spectra were 
taken from the ESO UVESPOP archive l|Bagnulo et al.l 
l2003h . The star HD 29907 was observed within the p roject 
67.D-0086A, as described by IMashonkina et al.l (l2003l l. De¬ 
tailed description of the spectra observ ed at the Lick obser¬ 
vatory is given bv ISitnova et al.l (l2015ll . 


4.2 Analysis of carbon atomic and molecular lines 

In this section we derive the carbon abundances of the se¬ 
lected stars from the atomic C I and molecular CH lines 
using the same approach as for Sun (Sect. [3]) and the same 
set of atomic and molecular data. 

The Shane/Hamilton observed spectra are affected 
by fringes in the IR spectral range. For the three stars, 
HD 59374, BD—4°3208, and BD-|-66°268, we applied the lo¬ 
cal continuum normalisation procedure and this made possi¬ 
ble to analyse the IR lines of C 1. For C I 9094 A and 9111 A 
in HD 103095 and HD 94028 we adopted equivalent widths 
from iTomkin et aP (ll992l L 

The NLTE calculations for C I were performed using 


0.3 
0.2 
X 0.1 
, 0 
U-0.1 
- 0.2 
-0.3 

-2.5 -2 -1.5 -1 -0.5 0 

[Fe/H] 

Figure 8. Differences in abundances derived from the atomic C I 
and molecular CH lines for the program stars. 



various S'h. As expected, everywhere NLTE leads to stronger 
C I lines and negative abundance corrections. For different 
stars Anlte varies between —0.06 dex and —0.45 dex for the 
IR lines and between 0.00 dex and —0.06 dex for the visible 
ones. The NLTE corrections decrease in absolute value to¬ 
wards lower effective temperature and metallicity. For exam¬ 
ple, for C I 9094 A in HD 59374 (5850/4.38/-1.02) Anltb 
= —0.34 dex (S'h = 0.3), and the corresponding value is 
—0.04 dex for the cooler star HD 103095 (5130/4.66/—1.26). 

The best fits of some atomic and molecular lines in 
HD 59374 and HD 84937 are shown in Fig.[71 

Carbon abundances obtained from the C I lines are pre¬ 
sented for all the stars in Table|6]and from the molecular CH 
lines in Tabled We comment below on the individual stars. 

For the most metal-rich star of our sample, HD 59374, 
both visible and near-IR lines of C I were measured, and we 
found the mean NLTE abundances from the two groups to 
be consistent, while Aloge(vis - IR) = —0.12 dex in LTE. 
The least NLTE differences, Aloge(vis - IR) = -1-0.07 dex 
and +0.05 dex, were obtained for S'h = 0.3 and 1.0, respec¬ 
tively. This is in line with our empirical estimate of Sh from 
analysis of the solar lines (Sect.[2|. When calculating the 
mean from the atomic and molecular lines in Table[7l we 
used everywhere the C I-based abundances determined with 
Sh = 0.3. 

For HD 94028 we analysed the single line in the visible 
region, C I 5052 A, and the two IR lines, with equivalent 
widths from iTomkin et ^ (Il992f l. 

Though the only atomic line, at 5052 A, was measured 
for HD 29907, the derived abundance is consistent with that 
from the molecular lines, independent of the line-formation 
scenario. 

For two IR lines of C I in HD 103095 we relied on 
equivalent widths from ITomkin et al.l (Il992ll and obtained 
an abundance difference of more than 0.2 dex between the 
lines, independent of the line-formation scenario. 

For two stars we cannot inspect the difference in abun¬ 
dances between the atomic C I and molecular CH lines. In 
BD—4°3208 the molecular CH lines are too weak to be mea¬ 
sured reliably, and only upper limit was estimated for the 
abundance from few individual lines (TablelT)). In contrast, 
in HD 122563 the atomic C I lines cannot be measured, 
but stellar carbon abundance is reliably obtained from the 
molecular CH lines, with a standard deviation of 0.04 dex. 

Figure[H] displays the abundance differences between the 
atomic C I and molecular CH lines for the Sun and seven 
metal-poor stars of our sample. Nowhere, Aloge(C I - CH) 
exceeds 0.07 dex and reveals any metallicity dependence. 
This suggests a minor influence of the 3D effects on abun¬ 
dance determinations from the molecular CH lines for metal- 
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Table 5. Characteristics of observed spectra. 


Object 

yi 

(mag) 

Telescope/ 

spectrograph 

Observing run 

Spectral range 

A 

R 

S/N 

HD 29907 

9.91 

8m VLT2/UVES 

Apr.2001 

3300-7000 

80 000 

150 

HD 59374 

8.50 

3m Shane/Hamilton 

Jan.2011, Mar. 2012 

3700-9300 

60 000 

100 

HD 84937 

8.28 

8m VLT2/UVES 

ESO UVESPOP 

3300-9900 

80 000 

200 

HD 94028 

8.22 

3m Shane/Hamilton 

Jan.2011, Mar.2012 

3700-9300 

60 000 

100 

HD 103095 

6.45 

3m Shane/Hamilton 

Jan.2011, Mar.2012 

3700-9300 

60 000 

100 

HD 122563 

6.20 

8m VLT2/UVES 

ESO UVESPOP 

3300-9900 

80 000 

200 

HD 140283 

7.21 

8m VLT2/UVES 

ESO UVESPOP 

3300-9900 

80 000 

200 

BD-4°3208 

9.99 

3m Shane/Hamilton 

Jan. 2011, Mar.2012 

3700-9300 

60 000 

100 

BD-|-66°268 

9.91 

3m Shane/Hamilton 

Jan. 2011, Mar.2012 

3700-9300 

60 000 

100 


^ V is a visual magnitude from the SIMBAD database. 







Figure 7. Best fits (continuum curves) of the selected atomic and molecular lines in HD 59374 (top row) and 
HD 84937 (bottom row). The observed spectra are shown by bold dots. The dashed curve in the right bottom panel 
shows the theoretical spectrum with no carbon in the atmosphere. 


poor stars. The only outlier is BD+66°268, with (Cl - CH) 
= 0.17 dex. 

We evaluated an influence of the uncertainties in stel¬ 
lar atmosphere parameters, namely 100 K in Tefj, 0.1 dex in 
log g, and 0.1 dex in [Fe/H], on the derived carbon abun¬ 
dances for five representative sets of Tefi/log g/[Fe/H]. The 
atomic lines were represented by C I 9094 A and the molecu¬ 
lar lines by the CH band 4310-4313 A. The obtained results 
are presented in Table[8l 


4.3 Carbon-enhanced stars 


Studies of the CEMP stars are important for better un¬ 
derstanding the Galaxy chemical evolution. High carbon 
abundance makes possible to measure lines of C I in the 
visible spectr al range even for v ery metal-poor (VMP) 
CEMP stars. Behara et al.l l|2010l . hereafter, BBLIO) and 
ISpite et al.l ll2013l . hereafter, SCB13) used C I 4932 A, 
5052 A, and 5380 A and the molecular CH lines to derive the 


carbon abundances of the four turn-off stars in the —3.3 ^ 
[Fe/H] ^ —2.5 metallicity range from the ID and 3D calcu¬ 
lations with applying the NLTE abundance corrections for 
lines of C I. They found large abundance discrepancies be¬ 
tween lines of C I and CH in both ID and 3D analysis, with 
(Cl — CH) of —0.35 dex to —0.79 dex in ID and —0.23 dex 
to -1-0.42 dex in 3D. 

We noticed that the NLTE carbon abundance correc¬ 
tion, Anltb = —0.45 dex, applied by BBLIO and SCB13 is 
rather large and performed the NLTE calculations using our 
model atom of C I and atmospheric parameters from BBLIO 
and SCB13. The obtained NLTE effects turned out much 
smaller, with Anlte = —0.04 dex. Using the ID-LTE abun¬ 
dances from BBLIO and SCB13 and our Anltes, we found 
consistent abundances from the two chemical species for 
three of four investigated stars, with (Cl — CH) = 0.02 dex 
(BBLlO’s star) and —0.03 dex and 0.07 dex (SCB13’s stars). 
For one of the BBLlO’s stars (Cl — CH) = —0.46 dex. In 
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Table 7. Carbon abundances of the program stars from the CH A-X electronic lines. 


HD/BD 

29907 59374 84937 94028 103095 122563 140283 -4°3208 +66°268 

A, A 

logec 


4362.4 - 4364.6 

6.77 

- 

6.58 

7.12 

7.07 

- 

6.35 

- 

6.34 

4366.0 - 4367.0 

6.78 

- 

6.58 

7.12 

7.10 

- 

6.32 

- 

6.39 

4310.0 - 4312.5 

6.73 

- 

6.60 

6.99 

- 

5.13 

6.35 

- 

6.34 

4313.4 - 4313.8 

- 

- 

6.64 

7.00 

- 

5.13 

6.30 

- 

6.39 

4218.710 

4218.734 

6.71 

7.48 

6.62 

7.08 

- 

5.09 

6.35 

- 

6.34 

4248.729 

4248.937 

4248.952 

6.74 

7.49 

- 

7.05 

6.93 

5.12 

6.38 

- 

6.34 

4253.000 

4253.206 

6.79 

7.55 

- 

- 

7.01 

5.09 

6.35 

- 

6.39 

4255.248 

4255.248 

6.74 

7.60 

6.67 

7.06 

6.93 

5.15 

6.40 

< 6.33 

6.34 

4263.969 

6.79 

7.57 

- 

7.05 

- 

5.01 

6.37 

- 

6.36 

4274.133 

4274.186 

6.75 

7.63 

- 

7.07 

- 

5.09 

6.42 

- 

6.36 

4356.355 

6.78 

7.64 

6.63 

- 

7.08 


6.46 

- 

6.37 

4356.371 

6.50 

- 


6.39 

- 

6.37 

4356.594 

6.76 

7.63 

- 

- 

7.05 

- 

- 

- 


Mean 

6.76 

7.57 

6.60 

7.06 

7.02 

5.10 

6.37 

< 6.33 

6.36 

a 

0.03 

0.06 

0.05 

0.04 

0.07 

0.04 

0.04 


0.02 

C I+CH 

6.78 

7.60 

6.60 

7.06 

7.01 

5.10 

6.40 

6.27 

6.45 

[C/Fe] 

-0.10 

0.19 

0.25 

0.13 

-0.16 

-0.71 

0.35 

0.04 

0.08 

Table 8. Changes in the abundances derived from the C I and CH lines caused by changes in stellar parameters. 

Models 

5100/4.5/-1 

5800/4.3/-1 

5800/3.7/-2 

6300/4.0/-2 4600/1.6/-2.5 


C I 

CH 

C I 

CH 

C I 

CH 

C I 

CH C I 

CH 

Areff=+100 K 

-0.17 

-0.01 

-0.08 

+0.12 

-0.06 

+0.18 

-0.05 

+0.17 -0.11 

+0.15 

Alog3=+0.1 

+0.04 

+0.01 

+0.03 

-0.01 

+0.04 

-0.03 

+0.03 

-0.04 +0.05 

-0.03 

A [M/H]=+0.1 

+0.03 

+0.10 

+0.03 

+0.06 

0.00 

+0.01 

0.00 

+0.01 +0.01 

+0.04 

Total 

0.18 

0.10 

0.09 

0.14 

0.07 

0.18 

0.06 

0.18 0.12 

0.16 


3D the difference (Cl — CH) is positive everywhere, up to 
+0.7 dex. 

This motivated us to compute the NLTE abundance 
corrections for lines of C I in the small grid of model atmo¬ 
spheres appropriate for the CEMP stars. We adopted Sh = 
0.3 and = 1.5 km s~ . The results are presented in Table[9l 
For the visible lines, C I 4932 A, 5052 A, 5380 A, and 6587 A, 
Anltb nowhere exceeds 0.13 dex in absolute value. The 
NLTE corrections are more negative, up to —0.70 dex, for 
the IR lines at 9061 A, 9078 A, and 9111 A. The SE calcula¬ 
tions were performed with [C/Fe] = 1, 2, a nd 3, although the 
empl oyed MARCS model atmospheres llGustafsson et all 
l2008h were computed with [C/Fe] = 0. To check how Anlte 


depends on the carbon abundance used in the model atmo¬ 
sphere constructio n, we asked F r ank Grupp to calcu late the 
two MAFAGS-OS (lGrupdl2004 iGrupp et al]l2009h models 
with common atmospheric parameters, 6250/4.0/—3, but 
different [C/Fe] = 0 and [C/Fe] = 3. When employing a 
common carbon abundance of [C/Fe] = 3 in the SE calcu¬ 
lations, we found very similar NLTE effects for C I in these 
two models. 

4.4 Comparison with previous studies 

iTomkin et al.l (1 19921 ') determined carbon abundances of 34 
metal-poor halo dwarfs from the atomic C I and molecular 
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Table 6. LTE and NLTE abundances of the program stars. 


HD, BD 

A, A 

LTE 

NLTE 




II 

p 

S'h^O.S 

29907 

5052.144 

6.82 

6.79 

6.79 

59374 

4932.049 

7.63 

7.59 

7.59 


5052.167 

7.75 

7.70 

7.71 


5380.337 

7.74 

7.70 

7.70 


Mean(vis) 

7.71 

7.66 

7.67 


cr 

0.07 

0.06 

0.07 


9062.492 

7.73 

7.51 

7.52 


9078.288 

7.80 

7.60 

7.62 


9088.515 

7.79 

7.59 

7.61 


9094.834 

7.88 

7.51 

7.54 


9111.809 

7.97 

7.68 

7.71 


Mean(IR) 

7.83 

7.58 

7.60 



0.09 

0.07 

0.08 


Mean(C I) 

7.79 

7.61 

7.63 


cr 

0.10 

0.08 

0.08 

84937 

8335.148 

6.69 

6.60 

6.62 


9062.492 

6.68 

6.59 

6.61 


9078.288 

6.68 

6.58 

6.60 


9111.809 

6.67 

6.55 

6.57 


Mean 

6.68 

6.58 

6.60 


cr 

0.01 

0.02 

0.02 

94028 

5052.167 

7.19 

7.14 

7.16 


9094 . 834 I 

7.18 

6.90 

6.94 


9 III. 8 O 9 I 

7.19 

7.01 

7.04 


Mean 

7.19 

7.02 

7.05 


cr 

0.01 

0.12 

0.11 

103095 

9094.8341 

6.92 

6.85 

6.88 


9111.8091 

7.14 

7.07 

7.10 


Mean 

7.03 

6.96 

6.99 


cr 

0.16 

0.16 

0.16 

140283 

8335.148 

6.43 

6.34 

6.35 


9061.433 

6.57 

6.47 

6.48 


9062.492 

6.52 

6.44 

6.45 


9078.288 

6.51 

6.41 

6.42 


9088.515 

6.55 

6.44 

6.45 


Mean 

6.52 

6.42 

6.43 


cr 

0.05 

0.05 

0.05 

-4° 3208 

9061.433 

6.34 

6.21 

6.22 


9062.492 

6.39 

6.29 

6.30 


9111.809 

6.44 

6.29 

6.30 


Mean 

6.39 

6.26 

6.27 


cr 

0.05 

0.05 

0.05 

+66°268 

9061.433 

6.64 

6.56 

6.58 


9062.492 

6.54 

6.47 

6.48 


Mean 

6.59 

6.52 

6.53 


cr 

0.07 

0.06 

0.07 


^ Equivalent widths from iTomkin et ^ lll992h . 


CH lines and found that the atomic lines lead to systemati¬ 
cally higher abundances compared with those from the CH 
lines, with the average difference, [C/Fe]ci— [C/Fe]cH = 
-1-0.4 dex. The C I lines were analysed based on the NLTE 
line formation. For four stars in common with that study 
we obtained [C/Fe]ci— [C/Fe]cH = —0.02 dex from our 
calculations and +0.25 dex from the data of iTomkin et al.l 


Table 9. NLTE abundance corrections for C I lines in the carbon- 
enhanced model atmospheres. 


A, A 

4932 

5052 5380 

6587 

9061 

9078 

9111 

log 9 


Teff = 4500 K, [M/H] = 

-3, [C/F( 

=] = 2 


0.5 

-0.10 

-0.12 -0.12 

+0.01 

-0.26 ■ 

-0.21 

-0.26 

1.0 

-0.08 

-0.09 -0.10 

0.00 

-0.21 

-0.17 

-0.21 

2.0 

-0.06 

-0.06 -0.07 

-0.02 

-0.13 ■ 

-0.10 

-0.13 



Tefj = 5000 K, [M/H] = 

-2, [C/F( 

,] = 2 


1.5 

-0.08 

-0.12 -0.10 

-0.06 

-0.55 ■ 

-0.50 

-0.59 

2.0 

-0.06 

-0.09 -0.07 

-0.06 

-0.48 

-0.41 

-0.50 

2.5 

-0.04 

-0.07 -0.05 

-0.05 

-0.37 

-0.32 

-0.39 



Teff = 5500 K, [M/H] = 

-2, [C/F( 

=] = 2 


2.0 

-0.08 

-0.13 -0.10 

-0.09 

-0.66 ■ 

-0.59 

-0.70 

3.0 

-0.05 

-0.06 -0.04 

-0.05 

-0.44 

-0.40 

-0.49 

4.0 

-0.02 

-0.03 -0.02 

-0.02 

-0.21 

-0.20 

-0.25 



Teff = 6250 K, [M/H] = 

-2, [C/R 

,] = 1 


3.0 

-0.02 

-0.02 -0.02 

-0.03 

-0.24 ■ 

-0.19 

-0.26 

4.0 

-0.01 

-0.01 -0.01 

-0.01 

-0.14 

-0.10 

-0.15 

5.0 

-0.01 

-0.01 -0.01 

-0.01 

-0.06 

-0.05 

-0.06 



Teff = 6250 K, [M/H] = 

-3, [C/R 

-] =3 


3.0 

-0.05 

-0.10 -0.07 

-0.07 

-0.65 

-0.55 

-0.69 

3.5 

-0.04 

-0.06 -0.05 

-0.05 

-0.54 - 

-0.45 

-0.56 

4.0 

-0.03 

-0.05 -0.03 

-0.04 

-0.41 

-0.34 

-0.39 

4.5 

-0.02 

-0.03 -0.02 

-0.02 

-0.29 ■ 

-0.24 

-0.29 

5.0 

-0.01 

-0.02 -0.01 

-0.01 

-0.19 

-0.14 

-0.19 


(I 1992 II . The discrepancy between the two studies can be due 
to using a 130 K, on average, lower effective temperature in 
iTomkin et al.l lll992h compared with our work. For the com¬ 
mon stars we employed th e IRFM effective temperatures, 
while iTomkin et al.l jl992ll the temperatures based on the 
b — y, R — Ij, R — Ik, and V — K colors. It can be seen from 
Table[8]that using the lower temperature leads to the higher 
C I-based and the lower CH-based abundance. 


5 CONCLUSIONS 

We constructed a comprehensive model atom for C I using 
the most up-to-date atomic data and evaluated the NLTE 
line formation for C I in classical ID models representing 
the atmospheres of late-type stars, where carbon abundance 
varies from solar value down to [C/H] = —3. 

In line with the previous studies we found that NLTE 
leads to stronger C I lines compared with their LTE 
strengths and negative NLTE abundance corrections. The 
deviations from LTE are large for the strong lines, C I 9061- 
9111 A (multiplet 3), 9405 A, and 9658 A, that form in 
the uppermost atmospheric layers, where collisions are in¬ 
efficient. For these lines Anltbs range between —0.10 dex 
and —0.45 dex depending on stellar parameters. The NLTE 
corrections do not exceed 0.03 dex in absolute value for the 
weaker visible lines C I 4932 A, 5052 A, 5380 A, and 6587 A. 
The NLTE effects depend strongly on the carbon abundance 
in the atmosphere. Therefore, the NLTE abundances of in- 
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dividual stars should be determined using an appropriate 
carbon abundance in the SE calculations. 


In the —1 ^ [M/H] 0 range our theoretica l re- 

sults agree we l l with that of IXakeda fc Hondal (|2005l i and 
iFabbiari et ahl (l2006l). At [M/H] = —2 our A nltes are con¬ 
sistent with that of iTakeda fc Hondal ll2005h . but they are 
0.08 dex smal l er in absolute value compared with that of 
iFabbian et al.l ll2006l i. We computed smaller NLTE abun¬ 
dance corrections for the [M/H] = —3 model atmo s phere s 
tha n the corresponding v alues in iTakeda &: Hond^ J2005l l 
and IFabbian et aP (l2006l i. This is not due to applying the 
iReidT i 1994 1 data. We did not find any observational data in 
the literature to check, what is correct. 


The treated model atom was applied to analyse lines 
of C I in the Sun and reference late-type stars covering the 
—2.56 ^ [Fe/H] $1 —1.02 range. The solar and stellar car¬ 
bon abundances were also determined fro m the molecular 
CH lin es, using recent molecular data from iMasseron et al.l 
(12014) . We aimed to check, whether different abundance in¬ 
dicators lead to consistent results for each star. 


The solar abundances derived from lines of C I, CH, and 
C 2 are consistent within the error bars, when applying S'h 
= 0.3 to 1 in the SE calculations for C I. Ou r results are in 
line w ith t he ID and 3D solar abundances of lAsplund et al] 
ll2009li an d ICaffau et al.l ll2010ll . 

For each program star the difference between the NLTE 
(Sh = 0.3) abundance from the atomic C I lines and the 
molecular CH lines does not exceed 0.07 dex. An exception 
is a cool VMP dwarf BD -1-66° 268, where the atomic IR lines 
are rather weak resulting in (C I - CH) = 0.17 dex. The 
obtained results suggest that the molecular CH lines can 
safely be used for stellar carbon abundance determinations 
with classical ID model atmospheres in a broad metallicity 
range, at least, down to [Fe/H] ~ —2.4. 

We show that the NLTE abundance corrections for 
C I 4932 A, 5052 A, and 5380 A in the carbon-enhanced 
6200/4.0/—3 models do not exceed 0.04 dex in abso¬ 
lute value. Applying our Anltb to the carbon ID-LTE 
abundances o f the four CEMP stars as determined by 
iBehara et al.l ll2010l ') and ISoite et al.l ll2013li leads to (C I — 
CH) of smaller than 0.07 dex for three of them. In 3D the dif¬ 
ference (Cl — CH) is positive everywhere, up to -|-0.7 dex. 
For the 4th star neither ID nor 3D helps to make abun¬ 
dances from the different chemical species consistent. We 
present the NLTE abundance corrections for lines of C I 
in the carbon-enhanced models of different Tefr, log g, and 
[M/Hj. They will be useful for accurate determinations of 
the carbon abundances of the CEMP stars. 

The NLTE method treated in this study for C I will be 
used in our further studies for accurate determinations of 
not only stellar carbon abundances, but also effective tem¬ 
peratures. The abundance difference between lines of C I 
and CH can serve as an efficient Teft indicator. 

Acknowledgements. We thank Frank Grupp for calculat¬ 
ing the MAFAGS-OS carbon-enhanced model atmospheres 
and Tatyana Ryabchikova for helpful discussion. This work 
was supported by the Russian Foundation for Basic Re¬ 
search (grants 14-02-31780 and 14-02-91153). 


REFERENCES 

Abate C., Pols O. R., Karakas A. I., Izzard R. G., 2015, 
A&A, 576, A118 

Akerman C. J., Carigi L., Nissen P. E., Pettini M., Asplund 

M. , 2004, A&A, 414, 931 

Alonso A., Arribas S., Martinez-Roger C., 1996, A&AS, 
117, 227 

Aoki W., Suda T., Beers T. C., Honda S., 2015, AJ, 149, 
39 

Asplund M., Grevesse N., Sauval A. J., Allende Prieto C., 
Blomme R., 2005, A&A, 431, 693 
Asplund M., Grevesse N., Sauval A. J., Scott P., 2009, 
ARA&A, 47, 481 

Bagnulo S., Jehin E., Ledoux C., Cabanac R., Melo G., 
Gilmozzi R., ESO Paranal Science Operations Team 2003, 
The Messenger, 114, 10 

Behara N. T., Bonifacio P., Ludwig H.-G., Sbordone L., 
Gonzalez Hernandez J. L, Caffau E., 2010, A&A, 513, 
A72 

Bensby T., Feltzing S., 2006, MNRAS, 367, 1181 
Brault J., Testerman L., , 1972, Preliminary Kitt Peak Pho¬ 
toelectric Atlas. Kitt Peak Nat. Obs., Tucson 
Brooke J. S. A., Bernath P. F., Schmidt T. W., Bacskay 
G. B., 2013, J. Quant. Spec. Radiat. Transf., 124, 11 
Burbidge E. M., Burbidge G. R., Fowler W. A., Hoyle F., 
1957, Reviews of Modern Physics, 29, 547 
Butler K., Giddings J., 1985, Newsletter on the analysis of 
astronomical spectra. No. 9, University of London 
Gaffau E., Ludwig H.-G., Bonifacio P., Faraggiana R., Stef¬ 
fen M., Freytag B., Kamp L, Ayres T. R., 2010, A&A, 514, 
A92 

Garbon D. F., Barbuy B., Kraft R. P., Friel E. D., Suntzeff 

N. B., 1987, PASP, 99, 335 

Casagrande L., Ramirez I., Melendez J., Bessell M., As¬ 
plund M., 2010, A&A, 512, A54 
Casagrande L., Schonrich R., Asplund M., Cassisi S., 
Ramirez L, Melendez J., Bensby T., Feltzing S., 2011, 
A&A, 530, A138 

Cohen J. G., Christlieb N., Thompson I., McWilliam A., 
Shectman S., Reimers D., Wisotzki L., Kirby E., 2013, 
ApJ, 778, 56 

Gollet R., Asplund M., Trampedach R., 2007, A&A, 469, 
687 

Greevey O. L., Thevenin F., Boyajian T. S., Kervella P., 
Chiavassa A., Bigot L., Merand A., Heiter U., Morel P., 
Pichon B., Me Alister H. A., et al. 2012, A&A, 545, A17 
Gunto W., Mendoza G., Ochsenbein F., Zeippen C. J., 
1993, Bulletin d’Information du Centre de Donnees Stel- 
laires, 42, 39 

Dobrovolskas V., , 2013, Effects of convection and non¬ 
equilibrium radiation transfer in stellar atmospheres. Vil- 
nus university. Center for physical sciences and technology 
Dobrovolskas V., Kucinskas A., Steffen M., Ludwig H.-G., 
Prakapavicius D., Klevas J., Caffau E., Bonifacio P., 2013, 
A&A, 559, A102 

Drawin H.-W., 1968, Zeitschrift fur Physik, 211, 404 
Fabbian D., Asplund M., Carlsson M., Kiselman D., 2006, 
A&A, 458, 899 

Froese Fischer C., 2006, Journal of Physics B Atomic 
Molecular Physics, 39, 2159 

Gonzalez Hernandez J. L, Bonifacio P., 2009, A&A, 497, 

















































14 S. A. Alexeeva & L. I. Mashonkina 


497 

Greggio L., 2005, A&A, 441, 1055 
Grupp F., 2004, A&A, 426, 309 

Grupp F., Kurucz R. L., Tan K., 2009, A&A, 503, 177 
Gustafsson B., Bell R. A., Eriksson K., Nordlund A., 1975, 
A&A, 42, 407 

Gustafsson B., Edvardsson B., Eriksson K., Jprgensen 

U. G., Nordlund A., Plez B., 2008, A&A, 486, 951 
Hayek W., Asplund M., Gollet R., Nordlund A., 2011, 

A&A, 529, A158 

Hibbert A., Biemont E., Godefroid M., Vaeck N., 1993, 
A&AS, 99, 179 

Huber K. P., Herzberg G., , 1979, Molecular spectra and 
molecular structure 

Ito H., Aoki W., Beers T. G., Tominaga N., Honda S., Car- 
ollo D., 2013, ApJ, 773, 33 

Jorgensen U. G., Larsson M., Iwamae A., Yu B., 1996, 
A&A, 315, 204 

Kennedy G. R., Stancliffe R. J., Kuehn G., Beers T. G., 
Kinman T. D., Placco V. M., Reggiani H., Rossi S., Lee 
Y. S., 2014, ApJ, 787, 6 

Kobayashi G., Karakas A. L, Umeda H., 2011, MNRAS, 
414, 3231 

KupkaF., Piskunov N., Ryabchikova T. A., Stempels H. G., 
Weiss W. W., 1999, A&AS, 138, 119 
Kurucz R. L., 1993, SYNTHE spectrum synthesis programs 
and line data 

Kurucz R. L., , 2007, Robert L. Kurucz on-line database 
of observed and predicted atomic transitions 
Kurucz R. L., Furenlid I., Brault J., Testerman L., 1984, 
Solar flux atlas from 296 to 1300 nm 
Lee Y. S., Beers T. C., Masseron T., Plez B., Rockosi C. M., 
Sobeck J., Yanny B., Lucatello S., Sivarani T., Placco 

V. M., Garollo D., 2013, AJ, 146, 132 

Luo D., Pradhan A. K., Saraph H. E., Storey P. J., Yu Y., 
1989, Journal of Physics B Atomic Molecular Physics, 22, 
389 

Lyubimkov L. S., Lambert D. L., Korotin S. A., 
Rachkovskaya T. M., Poklad D. B., 2015, MNRAS, 446, 
3447 

Mashonkina L., Gehren T., Shi J.-R., Korn A. J., Grupp 
F., 2011, A&A, 528, A87 

Mashonkina L., Gehren T., Travaglio G., Borkova T., 2003, 
A&A, 397, 275 

Masseron T., Plez B., Van Eck S., Golin R., Daoutidis 
L, Godefroid M., Coheur P.-F., Bernath P., Jorissen A., 
Christlieb N., 2014, A&A, 571, A47 
Moore C. E., 1972, A multiplet table of astrophysical in¬ 
terest - Pt.l: Table of multiplets - Pt.2: Finding list of all 
lines in the table of multiplets 
Otsuka M., Tajitsu A., 2013, ApJ, 778, 146 
Placco V. M., Frebel A., Beers T. C., Stancliffe R. J., 2014, 
ApJ, 797, 21 

Przybilla N., Butler K., Kudritzki R. P., 2001, A&A, 379, 
936 

Ralchenko Y. A., Kramida E., Reader J., Team N. A., 2008, 
NIST Atomic Spectra Database (version 3.1.5). USA 
Ramirez-Ruiz E., Trent! M., MacLeod M., Roberts L. F., 
Lee W. H., Saladino-Rosas M. L, 2015, ApJ, 802, L22 
Reddy B. E., Lambert D. L., Allende Prieto C., 2006, MN¬ 
RAS, 367, 1329 

Reetz J. K., 1991, Diploma Thesis. Universitat Miinchen 


Reid R. H. G., , 1994, Electron impact exitation of G I 
Rybicki G. B., Hummer D. G., 1991, A&A, 245, 171 
Seaton M. J., , 1962a, in Atomic and Molecular Processes 
(New York: Academic Press) 

Seaton M. J., 1962b, Proceedings of the Physical Society, 
79, 1105 

Shchukina N. G., Shcherbina T. G., 1990, Kinematika i 
Fizika Nebesnykh Tel, 6, 44 

Sitnova T., Zhao G., Mashonkina L., Ghen Y. Q., Liu F., 
Pakhomov Y., Tan K., Bolte M., Alexeeva S., Grupp F., 
Shi J.-R., Zhang H.-W., 2015, ArXiv:1506.01621 
Sluder A., Ritter J., Safranek-Shrader G., Milosavljevic M., 
Bromm V., 2015, ArXiv: 1505.07126 
Spite M., Gaffau E., Bonifacio P., Spite F., Ludwig H.-G., 
Plez B., Christlieb N., 2013, A&A, 552, A107 
Spite M., Cayrel R., Plez B., Hill V., Spite F., Depagne E., 
Frangois P., Bonifacio P., Barbuy B., Beers T., Andersen 

J. , Molaro P., Nordstrom B., Primas F., 2005, A&A, 430, 
655 

StarkenburgE., Shetrone M. D., McConnachie A. W., Venn 

K. A., 2014, MNRAS, 441, 1217 

Steenbock W., Holweger H., 1984, A&A, 130, 319 
Stuerenburg S., Holweger H., 1990, A&A, 237, 125 
Takeda Y., 1992, PASJ, 44, 649 
Takeda Y., Honda S., 2005, PASJ, 57, 65 
Tomkin J., Lemke M., Lambert D. L., Sneden C., 1992, AJ, 
104, 1568 

Urdahl R. S., Bao Y., Jackson W. M., 1991, Chemical 
Physics Letters, 178, 425 

Vernazza J. E., Avrett E. H., Loeser R., 1981, ApJS, 45, 
635 



